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Introduction 



Doppler surveys for exoplanets are monitoring almost every bright, chromospherically 
inactive solar type star (late F, G and early K types) and a fraction of the known M dwarfs 
within 30 parsecs. These successful programs have discovered about 400 exoplanets rang- 
ing in mass from a few times that of the Earth up to the lower mass limit for brown 
dwarfs (i.e., about 12 Mjup). This ensemble of exoplanets has revealed correlations be- 
tween both the metallicity and mass of the host star with the occurrance rate of exoplanets 
(ISantos. Israelian &: Mayoiil2004j : iFischer fc Valentill2005l : iLovis fc Mayorll2007l : I Johnson et al, 
2007h . 



The occurrence rate of g as giant planets ar ound M dwarfs appears to be lower than for 
F, G and early K type stars. lEndl et al.l (I20031) detected only one planetary companion in 



a sample of 100 M dwarf primaries. iJohnson et al.l ( 120071 ) carried out a statistical analysis 
and showed that late K and early M dwarf stars have a 1.8 ±1.0% occurrence rate for Jovian 
planets compared to 4.2 ±0.7% for solar-mass stars and 8.9 ±2.9% for higher mass subgiants. 

However, the protoplanetary disks of M dwarfs may still be a robust environment for the 
formation of lower mass planets. iForveille et al.l ( l2009l ) find that 30% of the twenty known 
planets ^ yith Msinz < O.l Mjnp orbit M dwarf stars. Consistent with this result, microlensing 
surveys (ISumi et al.ll2009l ) find far more planets with masses in the range between a few to 20 
Mffl. Thes e empi rical results are compatible with core accretion models of planet formation. 



Ida fc LinI (120051 ) predict that gas accretion will be quenched by gap formation in disks 
with smaller aspect ratios; the low surface density in the protoplanetary disks of M dwarfs 
implies a low f ormation proba b ility f or gas giants because massive cores don't form before 



gas depletion. iLaughlin et al.l (120041 ) also surmise that a significant population of "failed 



Jupiters" with cores of a few Earth masses should exist around M dwarfs. 

M dwarfs constitute the majority of stars populating our galaxy. Among the ~ 150 stars 
within 8 parsecs, about 120 are early-type M dwarfs, while only 15 are G dwarfs. Low mass 
stars are attractive targets for search es of terrestrial ma s s planets within th e habitable zone 
where liquid water might be present ( iTarter et al.ll2007l : iGaidos et al.l 120071 ). The habitable 
zone is closer to low luminosity stars and the shorter period planets induce larger, more easily 
detected reflex velocities. Yet, fewer than ten percent of l ate K and early M dwarf stars 



withi r i 30 parsec are b eing monitored by Doppler surveys (iMayor et al.l 120091 : iButler et al. 



20061 : lEndl et al.l l2003l ) . These stars are challenging targets because they are intrinsically 



faint at optical wavelengths. Indeed, many late type M dwarfs beyond ten parsecs are only 
now being detected; recently, the SUPERBLINK survey (lLepinell2005l : iLepine fc Sharall2005l ) 
identified a few thousand M dwarfs closer than 30 parsecs, opening up the possibility for 
larger scale surveys of these stars. 
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2. A New M and K dwarf Survey 

To learn more about the rate of exoplanet occurrence for M dwarf stars and mass 
and stellar metallicity dependences, we have started a survey of 1600 late K and early M 
dwarfs. The target stars are drawn from the SUPERBLINK survey, which provides an 
all-sky census of stars with proper motion larger than 40 mas yr~^ down to a magnitude 
limit R ~ 19.5. Main sequence stars are identified based on their location in a reduced 
proper motion diagram. A fraction of the stars are listed in the Hipparcos catalog; for 
those stars we calculate the distance from their Hipparcos parallax. All other stars have 
their d i stance s estimated based on the empirical V,V — J color-magnitude relationship from 



Lepind (120051 ). We select as probable late K and M dwarfs all stars with optical-to-infrared 
color V — J > 2.75, which roughly corresponds to absolute visual magnitudes My > 8.5 and 
masses M < 0.6Mq. We further narrow down our sample to stars with photometric distances 
d < 50 pc and visual magnitudes V < 12.5. The low proper motion limit and short distance 
range minimizes the kinematic bias, and the sample is expected to be complete for stars 
with transverse velocities vt > 9.5 km s~^ [vj- > 5.7 km s~^ for the 30 pc subsample), which 
effectively samples the local disk population over the full range of typical stellar ages, a s the 



selection even includes stars from nearby young moving groups iLepine fc SimonI (120091 ). Of 



the final sample of 1600 stars, 40% are listed in the Hipparcos catalog. 

We carry out low resolution spectroscopic screening of all candidates before they are 
observed at Keck using the Mark HI spectrograph at the 1.3-m McGraw-Hill Telescope at 
teh MDM Observatory at Kitt Peak. These observations are used to confirm the spectral 
type and main sequence status before beginning observations at Keck. 



Modeled after the N2K program (IFischer et al.l l2005l ) . our program uses a quick- look 



strategy to obtain three or four nearly consecutive observations of target stars to fiag short- 
period candidates. We then continue to obtain approximately log-spaced time series obser- 
vations to detect longer period planets. Here, we present the first planet to emerge from this 
program. 



3. HIP 79431 



3.1. Stellar Characteristics 



11. 34, J = 7.56, H 



HIP 79431 is an M3V star with apparent magnitudes of V 
6.86, K = 6.59 and color B — V= 1.486. The Hipparcos parallax (lESAlll997t Ivan Leeuwen 



20071 ) of 69.46 ± 3.12 mas yields a distance of 14.4 parsecs and we use this to calculate the 
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absolute magnitudes My 
fiRA = 28.31 mas yr~^ 



10.47, Mj = 6.69, Mh = 5.99, Mk = 5.72. The proper mo tions 
fJ-Dec = —208.36 mas yr~^ measured by the Hipparcos mission ( lESA 



19971 ) and our measured absolute radial velocity of -4.7 km s ^ correspond to the small space 
velocities oi U = +8.5 km s-\ V = +2.1 km s-\ W = -4.8 km s'^. 



Delfosse et al.l (120001 ) have established mass-luminosity calibrations for low mass stars 



and find low dispersion relations for Mj, Mh and Mx- The average of these three infrared 
mass- luminosity relations yield an average stellar mass of 0.49 Mq for HIP 79431. 

It is challenging to assess m etallicity froni M dw arf spectra, largely because of uncertain- 
ties in the molecular line data. iBonfils et al.l ( 120051 ) first der ived a photometr i c cali bration 
for M dwarf metallicities and this work has been updated by I Johnson fc Appj (120091 ). Both 
groups used spectral synthesis modeling to determine [Fe/H] for bright F and G type stars in 
binary systems with M dwarf companions and assigned the metallicity of the primary star to 
the M dwarf com panion. Based on t he sta r's height above the Mr- vs V — K main sequence, 
the calibration of I Johnson fc Appd (120091 ) yields a metallicity of [Fe/H] = +0.4 ±0.1. 



The effective temperature for M dwarfs is usually derived from a col o r-tern perature rela 
tion. The K-band temperature relation calibrated by ICasagrande et al.l (120081) is insensitive 



to me tallicity and M dwarf fluxes peak near the infrared K band. The ICasagrande et al. 
(l2008h relation has a scatter of only 19 K and yields Teff= 3191 K for HIP 79431. This is 



similar to the effective temperat u re of 3236 K that can be derived with the less robust B — V 
calibration of IValenti fc Fischerl ( 120051 ) for FGK type stars. 



Because the continuum flux near the Ca II H & K lines is so l ow in M dwarfs it is t ypica l 
to see strong emission from the core of the Ca II H & K lines. Ilsaacson fc Fischeii (120091 ) 
measure Shk = 1-15 ± 0.06 in this star. The logi?^K5 activity-based rotation period and 
ages a re not calibrated for stars with B — V greater than 1.0, however Ilsaacson fc Fischer 
(120091 ) find that this value of Shk is in the 50th percentile for stars of this color, so the star 
does not have an unusually high level of chromospheric activity for a dwarf star with this 
B — V color. The stellar characteristics are summarized in Table 1. 



3.2. Doppler Observations and Keplerian Fit 



HIP 79431 was added to the Keck program in April 2009 as part of our exoplanet survey 
of low mass stars. We obt ained 13 Doppler measurements of HIP 79431 over six months using 
the HIRES spectrometer (IVogt et al.lll994l ). Exposure times of 600 seconds for this y = 11.3 
star yielded a signal-to-noise ratio of just under 100. 
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Our Doppler analysis makes use of an iodine absorption cell in the light path before 
the entrance slit of the HIRES spectrometer. The iodine absorption lines in each program 
observation are used to model the wavelength scale a nd the instrumental profile of the 



telescope and spectrometer optics for each observation (IMarcy &: Butlerlll992l : iButler et al. 



19961 ). The typical Doppler precision for this faint star is about 3 m s~^. Based on the rms 
scatter of other stars on our program, we have empirically derived an additional error which 
may arise from astrophysical noise sources or systematic uncertainties in our data. Figure 
1 shows time series Doppler observations for stars with constant radial velocities and B — V 
colors similar to HIP 79431. The typical rms scatter is about 2.5 m s~^. We adopt 2.5 m s~^ 
as representative of the combined radial velocity jitter and systematic errors for late type 
stars and add this in quadrature with our formal single measurement errors when fitting and 
plotting our data. 

The observation dates, radial velocity data and measurement uncertainties for HIP 79431 
are listed in Table 2. In this Table, the measurement uncertainties are not increased by the 
2.5 m s~^ jitter term. The radial velocities exhibit an unambiguous signal. The radial 
velocities were fit with a Keplerian model and the uncertainties in the orbital parameters 
were derived using a bootstrap Monte Carlo analysis with 1000 realizations of the data 



(iMarcy et al.l l2005l ) in the following way: first, the theoretical velocities from the best-fit 
Keplerian model are subtracted, then the residual velocities are randomized, added back to 
the theoretical velocities and refit. This analysis captures the full range of uncertainties, in- 
cluding systematic errors, however, it will also scramble a signal from any additional planets, 
so it is only appropriate for good xl fit models. 

The Keplerian fit yields an orbital period of 111.7±0.7 days and a semi- velocity am- 
plitude of 149.5 ± 2.5 m s~^. The stellar mass of 0.49 Mq implies a companion mass with 
Msini = 2.1 Mjup and a semi-major axis of 0.35 AU. The parameters for the Keplerian 
orbit are listed in Table 3. 

The Doppler velocities are plotted in Figure 2 with a solid line indicating the best fit 
Keplerian model. In this plot, the single measurement errors have been increased by adding 
a modest 2.5 m s~^ for stellar jitter and systematic errors. The rms of the residuals to our 
Keplerian model is 3.9 m s~^ and the xl fit is 0.84, demonstrating that the fit is consistent 
with our single measurement errors and the added assumed jitter of 2.5 m s^^ As such a 
good fit implies, there is no evidence for additional planets. 
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4. Planet or Brown Dwarf? 

Because the Doppler technique derives M sin i and not the total mass, this leaves open 
the question of whether this is actually a stellar binary system in a nearly face-on orbit. 
HIP 79431 is only 14.4 pc away and a stellar companion on a nearly face-on orbit would 
induce motion around the barycenter of order 10 mas. The re-analysis of the Hipparcos 
astrometry is consistent (good ness-of-fit parameter of 0.15) with the prope r motion o f a 



single star ( Ivan LeeuwenI 120071 ). We used the intermediate astrometric data (lESAl 119971 ) to 
place an upper limit on the companion mass as a function of confidence level (one minus the 
false alarm probability FAP). The abscissa residuals to the Hipparcos parallax and proper 
motion solution span about 10 orbits and are plotted in Figure 3. There is no power in 
a Lomb-Scargle periodogram near the orbital period of 111 days in these data. Assuming 
gaussian statistics, a companion more massive than 0.3 Mq can be ruled out with 90% con- 
fidence (Figure 4), but, the Hipparcos data by themselves cannot eliminate the possibility of 
a late-type M or brown dwarf companion. We can, however, rule out the possibility that the 
location of HIP 79431 above the main sequence is due to a nearly equally-luminous compan- 
i on and therefore concl ude that this location is consistent with the high stellar metallicity 



(1 Johnson fc Appd 120091 ). 



Even a brown dwarf mass would require an improbable nearly face-on orbit. If we take 
the lower limit for stellar mass as 70 Mjup then the maximum inclination i for our system 
to surpass the stellar mass threshold is given by: 

Msini 

zro..p = arcsm,^^, (1) 

The probability that this particular inclination was observed can be calculated under 
the assumption that the orientation of orbits is random as: 



P{star) = cos(O) - cos(i7o jup) = 4 x 10 ^, (2) 

where i is the orbital inclination from Eqn (1). We can likewise calculate the probability 
that HIP 79431b falls in the brown dwarf range. This time, the inclination limit is far less 
severe since the planet only needs to have a true mass of 12 Mjup: 



Msini 



(3) 
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P{brown dwarf) = cos(i7o jup) - cos(ii2 jup) = 0.015 (4) 

So, in order for HIP 79431b to be a stellar companion, an inclination of less than 2° 
from face-on is required and that orientation has a probability of only 0.04%. If HIP 79431b 
is really a brown dwarf companion, the orbital inclination must be between 2° - 10°; the 
probability of this orientation is 1.5%. 

5. Transit Ephemeris 

Because HIP 79431b has a semi-major axis of 0.36 AU, the planet has a relatively low 
transit probability. The planet also resides in a moderately eccentric orbit that brings it as 
close as 0.25 AU to the host star at periastron. For the orbit of HIP 79431b, a; = 288 ± 3 
so a primary transit would occur when the planet is close to apastron. This further reduces 
the probability of a primary transit and we calculate a transit probability of only 0.5%. In 
this system, the secondary eclipse is somewhat more likely; without prior knowledge of a 
primary transit, we calculate a probability of 0.85% for the secondary eclipse. 

If a transit were to occur, the transit depth would be remarkable. A planet with a mass 
of 2.1 Mjup is supported by electron degeneracy and has a radius that is only a factor of four 
times smaller than the M dwarf host. We calculate a prospective primary transit depth of 
98 millimags and transit duration of almost 7 hours. The probability for a secondary eclipse 
is about double that for the primary transit, however the photometric decrement would only 
be about 0.004 millimags. 

Prospective transit times were predicted from the best-fit Keplerian parameters. The 
Thiele-Innes constants were used to project the orbit onto the plane of the sky and to 
determine the true anomaly corresponding to the central transit time, T^. The time of 
the next is calculated to be 2010 Feb 21 at UT 18:17:23:92 with an uncertainty of 1.6 
hours. The uncertainty in T^. was determined by fitting the 1000 trials of Keplerian fits (used 
to derive uncertainties in the orbital parameters). For each trial, the central transit time, 
ingress and egress were calculated and the standard deviation provided an estimate of the 
uncertainty. Future transit times can be estimated by rolling forward an integral number 
of orbital periods, however uncertainty in Tc increases with time from our last Doppler 
measurement. 
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Discussion 



As part of a new Doppler quick-look survey for late K and early M dwarfs, we have 
detected a planet with Msmi = 2.1 Mjup orbiting the M3V star, HIP 79431. The host star 
appears to be metal-rich with a photometrically-determined [Fe/H]= +0.4. 

Massive planets around late type stars are easy detections. However, the sample of 
M dwarf stars with planets is relatively small; only six have detected gas giant planets 
{Msini > O.SMjup). The remarkable system of planets orbiting GJ 876 (M4V) has three 
planets and two of these are gas giants in a 2:1 mean motion resonance; GJ 876b has Msinz 
= 1.935 Mjup and orbi tal period of about 30 d and GJ 876c has M si n i = 0.619 Mjup with an 
orbital period of 60 d jMarcv et al.lboOll . Il998l : lOelfosse et al.lll998h . GJ 849 (M3.5V) has a 
planet with Msinz= 0.82 Mjup in a 5-year orbit with significant residual velocities that now 
exhibit curvature, consistent with an additional planet in a wider orbit. GJ 317 (M3.5V) has 
a planet with Msini = 1.17 Mjup with a period of about 1.9 years and resi dual velocities 



consi stent with Jupiter-mass planet in an orbit of more than seve n vears (IJohnson et al. 



20071 ). GJ 832 (M1.5V) has a Jupiter-like planet in a 9.4-vear orbit JSailev. J. et al.l [2009[) 



GJ 179 (M3.5V) has a Jupiter-mass planet in a 6.3 ve ar orbit (IHoward et al.ll2009l ). GJ 649 
(M1.5V) has a Saturn-mass planet in a 1.6 year orbit ( Johnson et al. 20091 ). 



It is curious that half of th e known giant plane t s orb iting low mass M dwarf stars 



reside in relatively wide orbits. Kennedy fc KenyonI (120091 ) find that grain growth time 



scales are limited by photo evaporation with more rapid loss of disks in high mass stars. 
They note that the stellar mass-dependent disk dispersal timescales may account for the 
fact that planets orbiting higher-mass stars have wider orbits, since the disk may disperse 
more rapidly and halt inward migration. Perhaps lower surface density creates an analogous 
migration-hindering environment for low mass stars. Or, perhaps the protoplanetary disk of 
M stars forms fewer planets so that planet-planet scattering events are less common. With 
less competition for raw materials in the disk, it might be possible for single or well-spaced 
giant planets to dominate at wide separations. This is also a domain where gravitational 
instability could play a role boss06. 

HIP 79431b is one of the most massive planets detected in orbit around an M dwarf star. 
The calculated probability that this is actually an equal mass stellar companion is vanishingly 
small: 4 x 10~^. The The Hipparcos astrometric data do not rule out the possibility that 
this is a low mass star or brown dwarf companion, but we calculate only a 1.5% probability 
that HIP 79431b has an orbital inclination less than 10° and a true mass greater than 12 
Mjup. Statistically, this leaves a 98.5% probability that HIP 79431b is a planet with a mass 
between 2.1 and 12 Mj^p. 
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The Keplerian motion induced by a single planet can explain our data, although much 
less massive planets may exist on interior or exterior orbits. The excellent agreement of the 
single planet Keplerian planet solution with the data limit the mass of any second planet 
in the system. We performed a linear perturbation analysis around the one-planet solution 
and found that we can rule out planets more massive than M3 sin i of 12 or 15 Earths near 
the 3:1 or 2:1 interior resonances, respectively. 

Although Doppler surveys show that massive companions to M dwarf stars are rare, 
evidence for a few Jovian-mass planets candidates has emerged from microlensing surveys 
that sample wider orbits. (Because low mas s stars dominate th e galactic population, they 
are common lenses for background stars.) iGaudi et al.l ( 120021 ) place an upper limit and 
find that fewer than 33% of M dwarfs have Jupiter- mass companions between 1.5 and 4 AU. 
Neverth eless, only a small number of Jov i an- mass planets have been detected in microlensing 
surveys JSond et al.lbooi ludalskil [20051 ICaudi et al.lboosi : lOong et al.l[2009l ). 
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Fig. 1. — Doppler measurements of four M dwarf stars with constant radial velocities. These 
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Fig. 2. — Time series radial velocities from Keck Observatory are plotted for HIP 79431 
with 2.5 m s~^ of expected velocity jitter added in quadrature with the single measurement 
uncertainties. The Keplerian model is overplotted with an orbital period of 111 days, velocity 
amplitude of 149.5 m s~^ and eccentricity, e = 0.29. With these parameters and the stellar 
mass of 0.49 Mq, we derive a planet mass, Msinz = 2.1 Mjup and semi-major axis of 0.36 
AU. 
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Fig. 3. — Astromctric residuals (unbinned) from the Hipparcos intermediate data after fitting 
for parallax and proper motion. 
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Fig. 4. — False alarm probability (one minus the confidence level) of the companion exceeding 
a given mass given the Hipparcos abscissa residuals for HIP 79431 
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Table 1. Stellar Parameters: HIP 79431 



Parameter 


V 


11.34 


Mv 


10.47 


Mk 


5.72 


B-V 


1.486 


Spectral Type 


M3V 


Distance [pc] 


14.4 


Tcs 


3191 (1000 


[Fe/H] 


+0.4 (0.1) 


Mstar 


0.49 (0.05) 


'S'hk 


1.15 (0.06) 


Radial Velocity [km s~^] 


-4.7 



Table 2. Radial Velocities for HIP 79431 



RV (Trv 



JD-2440000 


(m s 




(m 


s-^) 


14928.07496 


-63 


.11 


2. 


,73 


14929.11611 


-70, 


.45 


2, 


,84 


14984.87477 


103, 


.92 


2. 


,68 


14986.96264 


125, 


.53 


3, 


,14 


14987.94866 


140, 


.47 


3, 


,10 


15028.96005 


-7, 


.06 


3. 


,88 


15041.89538 


-74, 


.88 


3. 


,21 


15042.92741 


-81, 


.75 


3. 


,55 


15048.75289 


-104, 


.55 


2. 


,80 


15074.72733 


-140, 


.47 


2. 


,84 


15078.74506 


-119, 


,78 


2. 


,90 


15084.73424 


-64, 


.40 


2, 


,97 


15106.73408 


139, 


.86 


3. 


,52 
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Table 3. Orbital Parameters for HIP 79431b 



Parameter 


Period (days) 


111.7 (0.7) 


Tp (JD) 


2454980.3 (1.2) 


u (deg) 


287.4 (3.2) 


Eccentricity 


0.29 (0.02) 


Ki (m s-^) 


149.5 (2.5) 


ttrel (AU) 


0.36 


Msini (Mj„p) 


2.1 


Nobs 


13 


RMS (m s-i) 


3.9 




0.84 



